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ABSTRACT 

A divergent geostrophic prediction model proposed by Duthie is 
investigated to determine its usefulness in producing realistic fore- 
casts and to determine the feasibility of further research using the 
gradient wind. 

24-hour forecasts were made using the prediction equation, 

۷ 25 - ۷ )لول + ( بوک رعال- - م۷ ۰ م22‎ 2 eof) 
(A) (B) HCA (D) 
at 300 mb, 500 mb, and 850 mb for summer and winter situations, with and 
without term B and with a tuning constant on term D. 

Tables presenting the error fields for these tests are included, 
along with the plotted forecast and verifying maps. These tables and 
illustrations will show that the best results at 500 mb and 300 mb are 
obtained with the B term out and a tuning constant of minus four and 
minus one respectively on the D term. At 850 mb different tuning con- 
stants produced no significant changes in the prognostic charts. 

The results of this investigation show that usable prognostic maps 
are produced, that the model tested has operational potentiality, and 


that further testing is warranted. 
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TABLE CF SYMBOLS 
the grid distance between consecutive points, 381km 
actual height of a surface minus the standard height 
the coriolis parameter 
the mean coriolis parameter, the value for 45N 


the upward component of the apparent gravitational 
acceleration 


the map factor 

millibar 

smoothed D values 

the twisting term, E: 22 
the horizontal vector wind 
the horizontal geostrophic vector wind 

the relative vorticity, ۰ ۷ x ۷ 

the geostrophic relative vorticity, ۴ ٧ x Vg 

the absolute vorticity, 7] z = + f 

the vertical wind component in the x,y,p,t coordinate system 
the horizontal Jacobian operator 

the finite-difference approximation of J(A,B) 

the finite-difference approximation of så 

the divergence of the wind 

the divergence of the geostrophic wind 

the horizontal Laplacian operator 

the finite-difference approximation of ve 


Fleet Numerical Weather Facility 


E Introduction. 

With the advent of the high speed digital computer, forecasting 
shifted from a highly subjective art accomplished by individuals using 
hand-drawn analyses and various prognostic methods, many of which had no 
theoretical foundation, to a more sophisticated method of solving the 
complicated differential equations governing the motíons of the atmos- 
phere. Many theoretical models have since been tested and the proficien- 
cy of the forecaster seems to have improved. However many of these 
models have experienced difficulty in moving troughs and systems with 
sufficient speed. 

The purpose of this paper is to introduce a divergent 5766 
model proposed by Duthie [2]. This model allows poleward moving parcels 
to increase their relative vorticity and equatorward moving parcels to 
decrease their relative vorticity, thus allowing for development and fill- 
ing. It will also move waves with a speed of the basic zonal current. 

The objectives are to show that this model, which uses the diver- 
gence of the geostrophic vind, will: 


a. numerically converge and produce a recognizable weather 


map; 
b. move waves at a speed consistent with the verifying analyses; 
and c. produce results sufficiently accurate to warrant further 
investigation. 


The author wishes to express his appreciation to Professor W. D. 
Duthie of the U. S. Naval Postgraduate School, on whose work this paper 


is based, for his guidance and encouragement in this investigation. In 
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addition, appreciation is also expressed to the personnel of the U. S. 
Fleet Numerical Weather Facility for their assistance in this project. 


Special thanks are extended to Mr. Leo Clarke for assistance in program- 


ming and many helpful suggestions. 


12 


2° Background. 


The isobaric vorticity equation for frictionless flow is [3] 
29 - - 7 ۷۰۱۷ + (250) 


It can be shown [1] that the divergence of the geostrophic wind 


is 


Vu = - Ng Sy رولا‎ f, (2.2) 


Assuming 


ad = 
GORE 


and substituting the geostrophic divergence for the divergence in equa- 


tion (2.1), the vorticity equation simplifies to 


Bet ) We (2.3) 


which expanded becomes 


pig E * 0 (2.4) 


Equation (2.4) shows that parcels moving poleward will increase 


their relative vorticity and vice versa. Expanding (2.4) further gives 
: a 
E =- gf L(=,39) * Dyd سل ج)‎ (1 (2.5) 


where 


33236 lv*z-ve: vef]. (2.6) 
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Substituting (2.6) into (2.5) gives the following 


Da عه‎ DE 
NS 0 vr ۰ ۷ ربا‎ - J (z, Sg) + وک‎ 0 (2) >, 


(A) (B) LE) (D) 


Equation (2.7) is a prognostic equation with a single variable, z. 
Term A is the local change, term B has the mathematical form of the 
Helmholtz term but with variable coefficients, term C is the advective 
term, and term D is the 9 term whose sign and magnitude now depend on 
the relative vorticity. 

It is proposed to test equation (2.7) with and without term B and 
with a tuning constant on term D. 

À basic linearized expression can be derived for equation (2.3). 
Assuming only a basic zonal current U and that products of perturba- 
tions are small and can be ignored, equation (2.3) becomes 


: 4 ۱ 
(2.8) 62 مدع 8 لله + 23 یہہ سے 


Now assuming a solution of the form 
7 ٥ 
equation (2.8) reduces to 
جو لاحو‎ 


This shows that waves will travel with the speed of the basic 


۱ 


current. This analysis shows that this model should move waves faster 


than models whose waves move with the speed of the Rossby waves [3], 
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3» Procedure. 
Data for 850, 300, and 500 mb were obtained from Fleet Numerical 
Weather Facility, Monterey, California, for summer and winter situations. 
These. data had been smoothed using a 15- pass, 5-point smoother over a 


3969-point grid (grid distance equal to 381 km at 60N) using equation (3.1). 
E 2 
S == DIN D, (3.1) 


where k is a constant smoothing factor. 


Using this smoothed field, equation (2.7) becomes 
a AS. 45 
Y dt VIS | Y Laf = -J (s, وک‎ ١+ ل نوك‎ Gs, bof), 
which after applying finite-difference approximations [5], becomes 
S | 
UL V*-(4»7- v).] RE == A 8 او ج) 0 ۱+24 جک ,دا‎ (3.2) 

where Sg = È ټپ ]| کسه‎ 85 - ۹ 5۰۲ ZF J ۰ 

Equation (3.2) was programmed for the CDC 1604 computer using a 
symbolic coded relocatable assembly program (SCRAP). 

The equation was solved for AS and then time stepped in increments 


of one hour using a forward time difference (3.3) the first hour and 


centered time differences (3.4) for succeeding hours. 


Str = Si FAS. (3.3) 
Seh  .و‎ AD, (3.4) 


At the end of 24 hours a prognostic map was plotted along with the 
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current analysis. From these two fields an error field was determined 
and from this error field a pillow (3.5) and RMSE (3.6)were computed for 


all points north of the equator. 


Pillow = 2 (A-B) (3.5) 


N 
RMSE = 


۶ ۱۸-۵۱-۲۳۱۶ GE 
N 


A represents the prognostic field, B the verifying field, and N the 


number of points used. 
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4. Results. 

Table 4.1 lists the pillow and the RMSE for tests run with the B 
term of equation (2.7) in and out for summer and winter situations. 

All runs were for 24 hours. This table shows that the B term tends to 
increase the error. 

Figure 4.1 is the prognostic map with the B term in. Superimposed 
on this map in dashed lines are the major troughs and centers of the 
prognostic map with the B term out (fig.4.2). Analysis of this map 
shows that the B term does not affect the movement of the waves but 
rather the intensity of the systems. The B term deepens lows and troughs 
and builds highs and ridges. The high over Northern Russia is 5349m 
with the B term in and 5321m with the B term out. It verifies as 5250m. 
The low over Greenland, Northern Canada, and the North Pacific covers 
a much broader area and is more intense with the B term in than the same 
low produced with the B term out. The verifying low is in better agree- 
ment with the low shown on the map with the B term out. 

Another matter which must be considered is the solution time on the 
computer. With the B term in, it takes 105 seconds to solve the equa- 
tion for one hour, 90 seconds of which is used to relax the "Helmholtz" 
term while it takes only 35 seconds to solve the equation with the B 
term out, 20 seconds of which is used in relaxation. Time T computer 
is a premium and must be considered in any model. 

Figure 4.2 is the prognostic map with the B term out. Superimposed 
on this map are the major troughs and centers of the verifying map 


(fig.4.3). Analysis shows that this model moves the waves faster than 
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B in vs B out 


500 mb 
Winter Summer 
Model P RMSE P RMSE 
B in -11' 266' 0! 102 
B out T 212" -2! 93! 





Table 4.1 RMSE Values for B term in vs. B term out, 500 mb 
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4.1 24-hour forecast for 500 mb from 8 December 1965, B term in 
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Figure 4,2 24-hour forecast for 500 mb from 8 December 1965, B term out 
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Figure 4.3 500-mb analysis, 9 December 1965 
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their verifying speeds. Centers also move too fast and are in general 
too intense. 

Recalling that in the perturbation equation (2.8), the e terms 
cancelled each other and thus resulted in the waves moving with a speed 
equal to the zonal wind, it was felt that by changing the sign of the e 
term in our model (D term) we would put some of the e effect back in 
and slow the waves down. 

Table 4.2 lists the error fields for summer and winter situations 
at 500 mb. This table shows that with a tuning constant of minus 4 on 
the D term, the model produces the best results. 

Figure 4.4 is the prognostic map with the B term out and a tuning 
constant of minus 4 on the D term. Superimposed in dashed lines are 
the major troughs and centers of the verifying map. 

Analysis of this map shows that the waves have slowed and now move 
with a speed consistent with the verifying waves. The centers appear 
to be in quite good agreement both in location and in intensity. It was 
55 that increasing the tuning constant resulted in continued slowing 
of the waves. 

It must be remembered that the verifying maps are prejudiced by 
FNWF's prognostic model since this is used əs a first guess to the verify- 
ing data. Therefore in sparse data areas (mostly open oceans) and wher- 
ever questionable data are received, the prognostic data are used. 

Figure 4.5 is a 500-mb prognostic map for summer with the B term 
out and a tuning constant of minus 4 on the D term. Figure 4.6 is the 
verifying map. Of interest in comparing these maps is the agreement on 


the trough coming out of the 5426m low east of Greenland and extending 
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B term out 


500 mb 
Winter Summer 
Tuning 8 Dec 11 Dec 3 July 10 July 
Constant P RMSE P RMSE P RMSE P RMSE 
1 L 212' 9' 2 -2' 93' 0! 80' 
-1 0, 160' 6' 176' -2' 87' -1' 76' 
-2 2۳ 13 a A 164 ' -2! 84! -1' 74! 
-4 3 107' -1' 160' : ۸ 81° 00 712 
-8 73 115' -4' 90! 30! 140! 


Table 4.2 Results of RMSE vs. tuning constant for 500 mb, B term out 
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Figure 4.4 24-hour forecast for 500 mb from 8 December 1965, B term out; 


tuning constant minus 4 
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Figure 4.5 24-hour forecast for 500 mb from 3 July 1966, B term out; 
tuning constant minus 4 
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Figure 4.6 500-mb analysis, 4 July 1966 
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down the European coast and out into the Atlantic west of Gibraltar. 
Again a trough in Western Russia which extends over the Black Sea and 
into the Eastern Mediterranean compares favorably. The low over Japan 
agrees in position but is 20m higher than the verifying height. The 
trough coming out of Alaska and down Western Canada into the United 
States does not show a closed low over North America as does the pro- 
gnostic map but a low of 5659m is indicated on the verifying map and 
could possibly have been drawn as a closed system. 

Table 4.3 lists the error fields for the 850-mb level for both 
summer and winter situations. This table shows a substantial decrease 
in the error values. Since this model is divergent, it might be expect- 
ed to perform better at a level where divergence has more effect. It 
should also be remembered that the gradients are not as strong as at 
500mb. 

At 850mb there does not seem to be any preference for a particular 
tuning constant. Figures 4.7, 4.8, and 4.9 are 24-hour forecasts made 
from 10 July 1966 with tuning constants of minus 1, minus 2, and minus 
4, respectively. Investigation of these maps will show no significant 
difference in location of the troughs and ridges or in intensity of the 
systems. Table 4.3 shows that varying the tuning constant does not ap- 
preciably change the RMSE value for any given day. 

Figure 4.10 shows a typical prognostic map for 850mb and figure 
4.11 is the verifying map. Investigation of these maps will show quite 
good agreement in both the location and intensity of the major systems. 

Test results, error fields, and investigations of the other maps 


in both summer and winter show similar patterns and substantiate these 
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B Term out 


850 mb 
Winter Summer 

Tuning 8 Dec 11 Dec 3 July 10 July 
Constant P RMSE P RMSE P RMSE P RMSE 
1 0 751 6' ۲ OR 48' 1 48! 
= 0 ' 69' 6' 68' وت‎ 47 ' 1 48 ' 
-2 0 67' 6' م٢‎ FR 48' qs 49 ' 
-4 p» 64' 6' 69! SA 49 ' p 49 ' 





Table 4.3 Results of RMSE vs. tuning constant for 850 mb, B term out 
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Figure 4.7 24-hour forecast for 850 mb from 10 July 1966, B term out; 
tuning constant minus 1l 
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4.8 24-hour forecast for 850 mb from 10 July 1966, B term out; 
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Figure 4.9 24-hour forecast for 850 mb from 10 July 1966, B term out; 


tuning constant minus 4 
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4.10  24-hour forecast for 850 mb from 8 December 1965, B term out; 


Figure 
tuning constant minus 4 


32 





سے س 


0 06۲ 09 DEC GS D 0دت‎ 850 MB HEIGHTS . 


۰ 


PROJECTION: POLAR STEREOGRAPHIC—TRUE AT 60’ NORTH LATITUDE FLEET NUMERICAL WEATHER FACILITY 


SCALE: 1:60,000,000 i DL . MONTEREY, CALIFORNIA © > 0 - 





Figure 4.11  850-mb analysis, 9 December 1965 
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initial results. 

Table 4.5 is a comparison between the error fields obtained with 
this model and the error fields of FNWF's model for the same dates at 
500 mb. The table shows that this model produces results comparable to 
those of the present operational model and in most cases might surpass 
it if this model were used as a first guess field for analysis. 

Since FNWF did not have an error field for 850mb, it was necessary 
to look elsewhere for comparison purposes. In 1962, Nicholson [4] in- 
vestigated a divergent model proposed by 'Arnason [1]. On the basis of 
the good results obtained from initial testing, 'Arnason's model was 
recommended for further investigation. Although an exact comparison is 
not possible since the days chosen for testing were not the same, it is 
felt that some benefit can be derived by comparing the order of magnitude 
of the error fields for each model. 

Table 4.6 shows the 850-mb RMSE values for both summer and winter 
situations. The verification scores for 'Arnason's model are those 
obtained using the optimum version of the model. The values entered for 
Duthie's model are those showing the largest errors. 

During initial testing at 300 mb, it was found that a relaxation 
coefficient of 17cm was to small. This coefficient required the computer 
to aake in excess of fifty passes to arrive at a solution. Since each 
pass takes approximately a second and a half, it required close to one | 
hour to make a 24-hour forecast including printing. It was decided to 
try 35cm on a trial basis. This cut the running time in half and was 
used throughout the remaining tests. However it is felt that this 
E 0 could Brill be increased and further experimentation is 


recommended. 
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Table 4.4 lists the RMSE values for tests run at 300 mb for both 
summer and winter situations. As was expected, it shows that the 
errors have increased substantially, being more than double those for 
500 mb. This can partially be explained by the large gradients found 
at this level. 

Figure 4.12 is the 24-hour prognostic map from 11 December 1965 with 
the B term out. Figure 4.13 is the verifying map. Investigation re- 
veals that the model is moving waves too rapidly eastward. Figure 4.14 
is the 24-hour prognostic map for the same time but with a tuning con- 
stant of minus one on the D term. Comparing figure 4.14 with figure 
4.12 shows that the tuning constant has slowed the waves down somewhat 
but they are still not in agreement with their verifying positions. 
Analysis of other maps with various tuning constants shows that none 
was able to compensate for the strong zonal wind at 300 mb. 6 
4.15 is the prognostic map for a summer situation with the B term out 
and a tuning constant of minus one on the D term. Figure 4.16 is the 
verifying map. In summer when the zonal wind is not quite as strong, 
the errors are reduced but the forecast waves are still moving too 
rapidly. 

It is felt that by using the gradient wind instead of the geostrophic 
wind in the C term of the model, it will reduce the speed of movement. 
Since the patterns of the prognostic and verifying maps are quite 
similar this modification should reduce the overall error field. 

During the course of this investigation it was found that the use of 
RMSE to define an error field did not do an adequate job. RMSE is only 


a measure of the difference in numerical values of two parameters. It 
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Summer 


3 July 1966 
RMSE 


223% 
216' , 
2227 


293 ' 


B term Out 


300 mb 


Winter 
Tuning 11 Dec 1965 
Constant P RMSE 
1 32 509! 
-1 1' 450' 
-2 -156' 5231 
-4 





Table 4.4 Results of RMSE vs. tuning constant for 300 mb, 


B term out 


Model vs Fleet Numerical 


500 mb 
Winter Summer 
8 Dec 11 Dec 3 July 10 July 
P RMSE P RMSE P RMSE P RMSE 
31 107' 1' 160' -3' 81' 0, 72' 
LI 105'18' 106' -9' 79' 2 79 


Model 


Duthie 
Model 


FNWF: 
Model 


Table 4.5 RMSE values for FNWF model vs. Duthie model, 500 mb 
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850 mb 


Winter Summer 
Model RMSE RMSE 
'Arnason 
"Modet. TII 87' 
Duthie 
Model 67' 50’ 





Table 4.6 RMSE values for 'Arnason model vs. Duthie 
Model, 850 mb 
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Figure 4.12 24-hour forecast for 300 mb from 11 December 1965 
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Figure 4.13 300-mb analysis, 12 December 1965 
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Figure 4.14 24-hour forecast for 300 mb from 11 December 1965, 
B term out; tuning constant minus 1 
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Figure 4.15 24-hour forecast for 300 mb from 3 July 1966, 


B term out; tuning constant minus 1 
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Figure 4.16 300-mb analysis, 4 July 1966 
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gives no indication as to how the maps resemble each other. 16 18 56 
that using RMSE in conjunction with a correlation coefficient, which 

is a measure of how much the products look alike, would provide a more 
useful error field. For example, a 300-mb map compared to a 500-mb map 
would show a high correlation coefficient but would also have a high 
RMSE. A high correlation coefficient together with a low RMSE would 
indicate that the prognostic map not only looked like the verifying 


map but also had small numerical errors. 
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55 Conclusions. 

On the basis of the results of this investigation, the model 
produces a prognostic map which is a very good representation of the 
verifying data. At both 500 mb and 850 mb the waves of the 'tuned' 
model move with a speed consistent with the verifying waves and the 
location and intensity of the centers agree remarkably well. At 300 
mb the model tends to move systems toofar. This is due to the use of 
a geostrophic wind with the strong zonal components at this level. 

In view of the close agreement at 500 mb and below, it is felt 
that further investigation is warranted, incorporating the gradient 
wind in lieu of the geostrophic wind, the geostrophic relative vorticity 
(equation 2.6 with the second term ín the brackets omitted) and daily 
testing begun over an extended period to determine the feasibility of 


this model for possible operational use. 
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